Direct emission of light mesons from quark-gluon plasma surface. 
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On the basis of hydrodynamic model of evolution of heavy ion collision we consider emission of 
7r, K, r), p, uj, K * mesons directly from the surface of quark-gluon plasma, with accounting of their 
absorption by surrounding plasma hadronic gas. We evaluate upper and lower limits on yields of 
these direct mesons in Pb+Pb collision at 158 A ■ GeV (SpS) and 3150 + 3150 A ■ GeV (LHC), and 
find, that even in the case of the lowest yield, direct r\ and heavier mesons dominate at LHC energy 
at soft pt (pt < 0.5 GeV/c). This enhancement of the low pt meson production can be considered as 
a quark-gluon plasma signature. 
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■ One of the goals of the ultrarelativistic heavy ion physics is a creation and registration of a new state of matter - 
quark-gluon plasma (QGP). However, being created in the heavy ion collision, QGP expands, cools and transforms 
back to hadronic gas, and up to now there is no reliable signal, which identifies uniquely its creation. 

As a rule, considering heavy ion collision with QGP creation, one assumes, that hadrons, produced on plasma 
surface, suffer numerous rescatterings in the surrounding hadronic gas, and thus final hadrons do not carry direct 
' information about plasma. In this paper we show that this widely accepted opinion is not quite right for the case of 
finite systems, created in the heavy ion collisions. Because of the comparability of free path lengths of hadrons in the 
hadronic gas with sizes of region, occupied by hot matter, significant number of hadrons, created on surface of quark- 
gluon plasma can pass through the surrounding hadronic gas without rescattering. Such kind of final state hadrons 
we call direct hadrons below. In contrast to usually considered freeze-out hadrons, originated due to evaporation from 
i the surface of hadronic gas or its freeze-out, direct hadrons carry immediate information about the plasma surface. 
£ — ' Emission of hadrons from the depth of hadronic gas was considered in literature several times. In the frame of 
quark-gluon string model it was shown that final hadrons, created in the heavy ion collision, originate not from 
thin freeze-out hypersurface, but from the whole volume, occupied by hot matter. Somewhat similar to our approach 
was used in papers ^] for evaluation of hadronic emission from the depth of hadronic gas. It was shown, that this 
' emission is not small enough to neglect, and can be one of the sources of enhancement of strangeness production 
Oh! in ultrarelativistic heavy ion collisions. Nevertheless, to our knowledge, until our previous paper [|3||, there was no 
special consideration of hadronic emission directly from the surface of QGP. We would like to stress, that in contrast 
to hadronic emission from the depth of hadronic gas, emission of direct hadrons from the surface of QGP gives unique 
' opportunity to explore QGP properties by means of hadrons. 
J> , In our previous paper || we considered S+Au collision at 200 A ■ GeV (SpS) and evaluated yields of pions, emitted 
copiously both from the surface of QGP and from the depth of hadronic gas before freeze-out. We compared their 
yields with yield of freeze-out pions and show that direct pions from QGP surface dominate in soft p t region, while 
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contribution of pions, emitted from the depth of hadronic gas is much smaller than contribution of freeze-out pions. 
Emission of direct pions from the surface of QGP result in enhancement of the yield of pions with low p t . This 
enhancement has been observed long ago in pt distributions of pions in AA collisions Q . A lot of explanations of this 
excess was proposed so far. Contributions of several proposed effects appear to be negligible small, but contributions 
of the strongest effects (decays of resonances ||, collective motion || and absence of chemical equilibrium Q) are 
sufficient for explanation of observed enhancement. Therefore, there is strong physical background for the direct pions 
from QGP. However, mentioned effects lead to much smaller contributions into distributions of heavier mesons. 

In this letter, in addition to pions, we consider emission of direct K, rj, p, oj and K* mesons in Pb+Pb collisions at 
158 A ■ GeV (SpS) and 3150+3150 A ■ GeV (LHC) energies and compare them with yields of corresponding freeze-out 
mesons. Estimation of situation at RHIC energies can be obtained by interpolation of our predictions at lower (SpS) 
and higher (LHC) energies. 

To estimate yields of the direct mesons in heavy ion collision we use the following model. Hot matter, created in 
the very beginning of collision, evolves hydrodynamically. On the background of this evolution the direct mesons 
are continuously emitted from the surface of QGP as a result of flying out of quarks and gluons from the depth of 
quark-gluon plasma, their hadronization on the plasma surface and fly out of the direct mesons through surrounding 
hadronic gas sometimes with rescattering. If direct meson suffers rescattering in the hadronic gas, then we assume, 
that it lost direct information about plasma and consider it further hydrodynamically. For freeze-out hadrons we 
assume thermodynamic and chemical equilibrium at freeze-out moment. 
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Probability for quark and gluon being emitted in the depth of QGP to reach its surface, and for meson to escape 
from hadronic gas without rescattering is determined by expression: 
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where integration is performed along (straight line) path of the particle in the hot matter, and A(e, x) - free path 
length of the particle, which depends on energy of the particle e and local energy density at the point x. We calculate 
free path lengths of quark and gluon in QGP and meson in hadronic gas using equation 
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where <xy(s) - total cross-section of interaction of i and j particles, mj and m 2 - mass of projectile and target 
particles correspondingly, T - temperature, and sum is taken over all possible two-particle reactions. Evaluating 
free path lengths of hadrons in hadronic gas we take into account only rescatterings on pions. In the case of the 
free path length of pions we use experimental cross-sections of mr scattering (see ||), while for K, 77, ... we assume 
a = OV+7T+ ~ 10 mbarn plus contributions from excitation of resonances. One could expect, that presence of nuclcons 
in the hadronic gas result in significant reduction of the free path length of the pion with respect to pure pionic 
gas (due to excitation of A resonances). However, we find, that for reasonable values of baryonic chemical potential 
fi < 200 MeV and temperatures below ~ 250 MeV, contribution of nucleons into free path length of pions is negligible. 

Emission rate (the number of particles, emitted from unit volume per unit time) of quarks and gluons from QGP 
we find from the condition, that infinitely thick layer of QGP emits quarks and gluons in accordance with Stephan- 
Boltsman formula. So we find: 
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where di - degeneracy, Xi - free path length of the particle, T - temperature. 

In this letter we are interesting in p t distributions at midrapidity region, therefore we describe evolution of the hot 
matter by use of Bjorken hydrodynamics with transverse expansion. To investigate sensitivity of our predictions to 
the description of freeze-out process, we consider two extreme versions of hydrodynamic model with pressure outside 
freeze-out surface P ou tside = and P ou tside = Pinside- 

As far as it is not possible to describe consistently hadronization of quarks and gluons on the plasma surface, we 
estimate upper and lower limits on the yields of direct hadrons, by use of two models of hadronization - 'creation' 
model, giving the lower limit and 'pull in' model, used in the paper || and giving the upper limit. In both models 
we assume, that a quark or gluon, flown through the plasma surface, pulls tube (string) of color field. Creation of the 
final hadron corresponds to the breaking of the tube due to discoloring of the moving out quark or gluon. In the first 
model we assume, that this discoloring takes place as a result of creation of a quark-antiquark pair in the strong field 
of the tube. This assumption is used in the well known family of event generators JETSET ||, where all parameters 
are chosen to fit the e + e _ annihilation at y/s = 30GeV. Therefore, in the 'creation' model we extract corresponding 
dependence from JETSET 7.4. In the second model we assume that discoloring takes place as a result of 'pulling in' 
of soft quark or gluon with corresponding color from the pre-surface layer of QGP into the tube. Because of the large 
number of soft quarks and gluons in the QGP each moving out quark and gluon can transform into some hadron if 
its energy is larger than mass of this hadron. If several hadrons can be formed, then we take their relative yields from 
JETSET 7.4. Having probabilities of hadronization of quarks and gluons R^ 9 (e q ), evaluated using these two models, 
we obtain the following hadronization function: 
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where s, 9, 4> energy, polar angle and azimuthal angle of initial quark (q) or gluon (g) or final hadron (h) correspond- 
ingly, and rrih - mass of the hadron. This fragmentation function is normalized to describe fragmentation of one quark 
or gluon to R'^ 9 {e q ) hadrons with energy in the range rrih < £ n < e q - 

To apply our model for central Pb + Pb collisions at 158 A ■ GeV we choose initial conditions (initial energy 
density e°„, time of thermalization Ti n and initial radius Ri n ), QGP-hadronic gas transition temperature (T c ) and 
freeze-out temperature (T/ reeze ) to reproduce experimental pt distribution of 7r° at midrapidity ||. Within two 
versions of the hydrodynamic model, mentioned above we use: = 2.5 GeV j fm 3 , Tj n = 1/m/c, Ri n = 7 fm, T c — 
170 MeV, T freeze = 130 MeV for P outslde = P msld e, and e° ln = 0.35GeV/fm 3 , r rn = 4.5 fm/c, R m = 10.5 fm, T c = 
150 MeV, Tf reeze = 140 MeV for P ou tside = 0. As far as in the case of P ou tside — the significant pressure gradient 
on the freeze-out surface result in strong acceleration of freeze-out hadrons, we choose lower initial energy density and 
transition temperature and higher freeze-out temperature than in the case of Poutside = Pinside- 

Evaluated pt distributions of various mesons at midrapidity for the 'creation' (left plot) and 'pull in' (right plot) 
models of hadronization of quarks and gluons for hydrodynamics with P ou tside — Pinside are shown on the fig. 1. For 
convenience, distributions of K, rj, ... mesons are subsequently multiplied by 10 -2 , 10 -4 , .... Solid lines correspond to 
direct mesons from QGP surface, and dashed lines - to freeze-out mesons. We do not distinguish ir + , ir°, and 7r~, 
K + , K" K°, and K°, so that we show distributions, averaged over all pions, kaons etc. Probability of fragmentation 
of low energy quarks and gluons in the 'pull in' model is higher than in the 'creation' model, therefore the number 
of direct mesons is higher on the right plot. In the both models of hadronization contributions of direct mesons 
increase with increasing of the meson mass, and in the case of rj and heavier mesons, direct mesons dominate over 
freeze-out ones, for the 'pull in' model. Unexpected, that direct mesons, emitted from the hottest region of the 
collision, contribute mainly into soft pt region. This takes place because during hadronization a quark or gluon 
'decays' onto direct hadron and part of the tube which pulls back to the plasma, so that the temperature of direct 
hadron significantly decrease with respect to typical temperatures in the system. 

The same distributions, evaluated in the case of Poutside = are shown on the fig. 2. In this case because of the 
faster fly out of the hadronic gas the width of the layer of hadronic gas, surrounding QGP, becomes smaller and it 
less screens direct hadrons. As a result, the contributions of direct hadrons become more significant: direct K, rj and 
heavier mesons dominate over freeze-out mesons for both models of hadronization. 

To investigate dependence of the yields of direct mesons on the energy of collision, we evaluated these yields in 
the central Pb+Pb collisions at LHC energy. We consider case Poutside = Pinside and use the same values of T c and 
Tfreeze as for SpS energies and the following initial conditions: e\ n = 20 GeV j fm 3 , r in — 1 fm/c, R in = 7 fm, what 
corresponds to the multiplicity at midrapidity dN* /dy ~ 2 • 10 3 , predicted by cascade models. Distributions of direct 
and freeze-out hadrons, evaluated in this case are shown on fig. 3. In contrast to SpS at LHC energy direct mesons 
dominate in the soft pt region for K and heavier mesons for both models of hadronization. 

As one can see, contribution of direct mesons enhances the pt distributions in the soft pt region. This takes place 
because during hadronization of quarks and gluons on the surface of QGP a part of their energy is pulled back to the 
QGP and thus direct hadrons have lower temperature then typical temperature in the system. As far as contribution of 
direct hadrons increase with increasing of their mass, we predict increasing of the low p t enhancement with increasing 
of the hadron mass - in contrast to all other possible explanations, such as decay of resonances, collective motion, 
absence of chemical equilibrium etc. This behavior can be used as signature of QGP creation. 

To conclude, we consider light mesons, emitted directly from the surface of quark-gluon plasma and escaped from 
hot matter without rescattering. This mesons gives unique opportunity to test the pre-surface layer of QGP via 
strongly interacting particles. We show, that direct mesons contribute mainly into soft pt region (pt < 0.5 GeV/c). 
Using two extreme models of hadronization of quarks and gluons on the plasma surface and two extreme models of 
hydrodynamic expansion, we estimate lower and upper limits on the yields of direct mesons with different masses. We 
find, that at SpS energies in the frame of both models of hadronization direct mesons dominate in the soft pt region 
for rj and heavier mesons in the case of Poutside = but only within the 'pull in' model of hadronization in the case 
of Poutside = Pinside- However, for LHC energy direct rj and heavier mesons dominate at soft p t for all considered 
models. Contribution of direct mesons result in appearing of the enhancement of pt distribution in the soft region, 
which can be used as a new signature of the of quark-gluon plasma creation. 
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FIG. 1. Yields of direct (solid lines) and freeze-out (dashed lines) mesons in Pb+Pb collision at SpS energy within 'creation' 
(left plot) and 'pull in' (right plot) models of hadronization for hydrodynamics with P outS ide = Pinside- Distributions of K, r), ... 
are multiplied by 10~ 2 , 10~ 4 , .... 



FIG. 2. The same as fig. 2, but for hydrodynamics with P out side = 0. 



FIG. 3. The same as fig.l, but for Pb+Pb collision at LHC energy. 
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